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Abstract 

Grover's search algorithm can be applied to a wide range of problems; even problems not 
generally regarded as searching problems, can be reformulated to take advantage of quantum 
parallelism and entanglement, and lead to algorithms which show a square root speedup over 
their classical counterparts. 

In this paper, we discuss a systematic way to formulate such problems and give as an example 
a quantum scheduling algorithm for an R\\C ma x problem. R\\C m ax is representative for a class 
of scheduling problems whose goal is to find a schedule with the shortest completion time in an 
unrelated parallel machine environment. 

Given a deadline, or a range of deadlines, the algorithm presented in this paper allows us to 
determine if a solution to an R\\C ma x problem with N jobs and M machines exists, and if so, it 
provides the schedule. The time complexity of the quantum scheduling algorithm is 0(V M N ) 
while the complexity of its classical counterpart is 0(M N ). 



1 Introduction 

In recent years we have witnessed significant theoretical and some encouraging experimental results 
in the area of quantum computing. In 1994, Peter Shor found a polynomial time algorithm for 
the factorization of n-bit numbers on quantum computers |2()| . His discovery generated a wave of 
enthusiasm for quantum computing, for two major reasons: the intrinsic intellectual beauty of the 
algorithm and the fact that efficient integer factorization is a very important practical problem. 
The security of widely used cryptographic protocols is based upon the conjectured difficulty of the 
factorization of large integers. 

Shor's algorithm reduces the factorization problem to the problem of finding the period of a 
function, but uses quantum parallelism to find a superposition of all values of the function in one 
step. Then the algorithm calculates the Quantum Fourier Transform of the function, which sets the 
amplitudes into multiples of the fundamental frequency, the reciprocal of the period. To factor an 
integer, Shor's algorithm measures the period of the function 1 . 

In 1996, Grover described a quantum algorithm for searching an unsorted database containing 
N items in a time of order y/~N while on a classical computer the search requires a time of order N 
|H] . The critical aspect of a quantum algorithm is to create a superposition of all possible states and 
amplify the "solution" . The speedup of Grover's algorithm is achieved by exploiting both quantum 
parallelism and the fact that in quantum theory a probability is the square of an amplitude. Bennett 
and his co-workers pQ and Zalka [2] showed that Grover's algorithm is optimal. No classical or 
quantum algorithm can solve this problem faster than time of order yN. 

A powerful version of the technique used by Shor is the phase- estimation algorithm of Kitaev 1121 



Grover's search can be applied to a large number of unstructured problems and lead to a square 
root speedup over the corresponding classical algorithms. For well structured problems classical 
algorithms that produce an approximate solution and perform faster exit. 

Recently, Furrow discussed applications based upon Grover-type algorithms [Sj. He considered 
three classes of applications: a) Graph algorithms, e.g. Breadth-First Search (BFS), Depth-First 
Search (DFS), bipartite matching, b) Computational Geometry algorithms, e.g. Maximum points on 
a line, c) Dynamic programming algorithms, such as coin changer. The author reports the quantum 
versus classical complexity for BFS and DFS, 0(\/VE lg V), versus O(E), with V the number of 
vertices and E the number of edges of the graph; for bipartite matching, 0(Vy/(E + V) lg V), versus 
0((E + V)y/V))\ for maximum points that lie on a line in R 2 , out of N points, 0(N 3 / 2 IgN) versus 
0{N 2 \gN), and so on. 

Most of the problems discussed in [S] are intrinsically search problems and the idea of applying 
Grover's search comes naturally to mind. There is an even larger class of problems which, at the first 
sight, do not seem directly related to Grover's search. Applications such as scheduling and resource 
allocation are not naturally search problems; nevertheless they share some common properties, can 
be reformulated to take advantage of quantum parallelism and entanglement, and lead to algorithms 
which show polynomial speedups over their classical counterparts. 

A scheduling problem is characterized by a tuple (a | (3 | 7) where a denotes the machine environ- 
ment, (3 summarizes the set of constraints, and 7 denotes the optimality criterion. The makespan 
of a schedule, C max is the maximum completion time of any job in the schedule. For example, 
-P||C maa: and R\\C max require the shortest makespan and apply to identical machine environment 
and, respectively, a non-homogeneous one. 

When we turn our attention to problems when a deadline is imposed, or when we wish to find 
a schedule with a given range of possible average completion time we discover that a full range of 
scheduling problems have a quantum counterpart which can take advantage of Grover's search. 

We illustrate these ideas with an example: given a deadline, or a range of deadlines, the algorithm 
presented in this paper allows us to determine if a solution to an R\ \C rnax problem with N jobs and M 
machines exists, and if so, it provides the schedule. The time complexity of the quantum scheduling 
algorithm is 0(V M N ) while the complexity of its classical counterpart is Q(M N ). 

Real-time systems are subject to deadlines and Quality of Service (QoS) constraints imposed to 
many systems require a given range of average completion times. Thus, the classes of scheduling 
algorithms we discuss in this paper are of significant practical importance. Such algorithms have a 
quantum counterpart that enjoy a square root speedup. 

2 Scheduling Algorithms 

Scheduling is the problem of assigning tasks to a set of resources subject to a set of constraints, over 
time, in an "optimal" manner. 

We are given a set of N = 2™ jobs, J — {J\,J%, ■ ■ -Jn} and a set of M = 2 m machines, 
M. = {M.\, M.2, ■ ■ ■ M-m}- A schedule S for the sets (J,M) specifies which units of time 
machine A4j uses to process job Ji. We call Cf the completion time of job Ji under the schedule 
iS. The makespan of the schedule iS is the maximum completion time of any job in schedule S: 

Cmax — max • 

We are often interested in scheduling problems involving multiple machines. We distinguish three 
cases in the parallel machine environment: 

1. Identical parallel environment. All the machines are identical and job Ji requires Tj units of 
time on any machine, 

2. Uniformly related parallel environment. Each machine M.j has a speed Sj > and job Ji if 
processed entirely on machine Mj would take Ti/sj units of time, and 
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3. Unrelated parallel machine environment. Different machines have different capabilities and the 
speed of machine Mj on job Ji is sy. Then the processing time of job Ji on machine Mj is 
Tij = Ti/sij, < Tij < Q and Q = V. 

These machine environment are denoted by P, Q, and R, respectively. 

Examples of scheduling constraints include deadlines (e.g., job i must be completed by time t), 
resource capacities (e.g., there are only 5 machines), precedence constraints on the order of tasks 
(e.g., one must be done before another), and priorities on tasks (e.g., finish job j as soon as possible 
while meeting the other deadlines). A priority rule assigns to job Ji a priority m. A busy schedule 
defines the situation where when one machine becomes available it starts processing the job with 
the highest priority. 

Each scheduling problem could have problem specific optimization criteria. The one which re- 
quires the minimization of the makcspan is referred to in scheduling theory as C m ax- Other opti- 
mization criteria can be considered. For example, we may wish to optimize the average completion 
time of all jobs: 

1 N 

i=l 

an optimization criterion denoted as X^iLi Cf . 

Many scheduling problems are MV— hard on classical computers [3| and have proved to be very 
difficult even for a relatively small instance. For example, a 10-job-10-machine job-shop scheduling 
problem posed in 1963 remained unsolved until 1989 [1J. Most classical scheduling algorithms gain 
speedup only on some special cases with relatively small instance or they try to find good schedules 
instead of the optimal one. Furrow|Sj also gave a special case of the P\\C ma x problem that can be 
solved by dynamic programming, which requires that the number of different jobs processing times 
be bounded by a constant. It turns out that there are polynomial approximations for some P\\C m ax 
problems. 

We consider an R\\C max scheduling problem in which all machines are unrelated. All jobs are 
available at the beginning time and there are no precedence constraints. As no preemption is allowed, 
once job Ji started processing on machine Mj it must complete its execution before another job, 
Jk can be processed on Mj. 

Given a set of N = 2 n jobs and M = 2 m machines we can construct 2 m2 ™ different sched- 
ules. This R\\C max scheduling problem is MV— hard and is difficult for classical algorithms even 
with small instance. Up to now, most classical algorithms use linear programming based rounding 
techniques (LP-rounding) to search an approximate schedule instead of the optimal one [81 119j. If 
the number of machines m is part of the input, the best approximation algorithm to date is a 2- 
approximation by Lenstra, Shmoys and Tardos |13j which can find a schedule with C max < 2* C^ x . 
Moreover, the problem cannot be approximated within a factor strictly smaller than 3/2, unless 
P=NP J3| . No proper classical algorithm addresses a general case of searching the optimal schedule 
except the exhausting search which has time complexity 0{M ) on a classical computer. The paper 
suggests a reformulation of such problems to take advantage of Grover-type search and gain square 
root speedup on searching the optimal schedules over their classical counterparts. 



3 Information Encoding 

Consider a set of N = 2" jobs running on M — 2 m unrelated machines. We assume that the jobs 
could have different processing times on different machines and that the processing times are integers 
in the range < T i3 < Q = 2 q , 1 < i < 2", 1 < j < 2 m . We also assume that we have a quantum 
system with r — m + q qubits for each job. 
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Given job Ji running on machine Mj, we encode the job-machine information as a vector | e\) 
obtained as the tensor product of the machine index, j, and the processing time, T^ : 

I 4) =1 J>® I T a)- 

Then we define job state vectors for job i, \ Ji) as any superposition of its job-machine vectors. First, 
we consider an equal superposition of job-machine vectors as: 

A schedule is a tensor product of job-machine vectors: 

\S k )=®Zi\4), 1<J,<2™, 
which includes one job-machine vector for each job. A specific machine may be present in multiple 
job-machine vectors, when the schedule requires multiple jobs to be executed on the same machine, 
or may not appear in any job machine vectors of the schedule | <Sfe) if none of the jobs is scheduled 
on that machine. 

The equal superposition of all schedules is: 



v k=i 



Let fi be an operator which given a schedule constructs the running time on each machine for 
that schedule. When applied to the equal superposition of all schedules, it produces a superposition 
of the running time | T) on each machine for all schedules: 

I ST) = fi(|S)|0)). 

where, | ST ) denotes the entangled state of S and T, while the tensor product of S and T is | ST) . 
Let A be an operator which computes a superposition of the makespan of all schedules: 

= A(| 5T) | 0)) = A 0(| 5) |0> |0)). 

Figure H outlines our procedure to produce an equal superposition of the makespans of all sched- 
ules. We now turn our attention to the quantum circuits to carry out the transformations discussed 
in this section and to the question how to obtain from the superposition of all makespans the optimal 
one. 

First, we need to prepare the job vector | J,-) in an equal superposition state which includes the 
processing times of job Ji on all machines, as shown in Figure |21 We use m index qubits to control 
the job-machine information encoding. As each index qubit is prepared in state 1 / v2(| 0)+ | 1)), 
the target qubits will be prepared in superpositions of all possible job-machine states, e\, 1 < i < 
n, 1 < j < m. 

Example: Tablensummarizes the processing time of 8 jobs on 4 machines, where < < 2 4 = 16. 
Thus n = 3, m = 2, and q = 4. The running time of J\ on machines Mi, M.2, M.3 and Mi are 
respectively, 1, 3, 7, and 15 units of time. 

The four vectors used to encode the processing time of job Ji on machines A^i, A^2, A^3 and 
Ma, are respectively: 

I ei) =| 000001), I ef) =| 010011), | e?) =| 100111), and | e\) =\ 111111). 
For J 2 the basis vectors are | 000010), | 010001), | 101001), | 110011), and so on. 
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Figure 1: Outline of our algorithm to prepare the makespan vector. First we prepare job vectors 
in an equal superposition of all possible schedules | S < -- N ^ n+q ^). Using a quantum circuit f2 with 
inputs | S' N ( m+q '') and a register of M(n + q) qubits in state | 0), we construct the superposition 
of the running time on each machine for every schedule. Then we construct the makespan of each 
schedule using the operation A. Superscripts indicate the number of qubits for vectors. Note that 
the number of jobs is N = 2", the number of machines is M = 2 m , and the maximum execution 
time of a job on any machine is Q = 2 q . 



Table 1: The processing time of 8 jobs on 4 machines 
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Figure shows the circuit used to prepare the job vector | J\) for our example. The job state 
| Ji) is prepared in an equal superposition of basis states: 

I Ji) = \{\ 000001)+ 1 010011)+ 1 100111)+ 1 mm)). 

We can prepare other job vectors in the same way, for example: 

I Ji) = g(| 000010)+ | 010001)+ | 101001)+ | 110011)). 

A schedule vector | S) is the tensor product of all job vectors. As each job vector is prepared as 
an equal superposition, the schedule vector is in the equal superposition of all possible schedules. 
We now provide two examples of schedule vectors: 

(i) Schedule S%: 



[Jx^Mi, J 2 ^M 2 , Js^Mi, Ji^M 4 , J 5 ^M 3 , J 6 ^M 2 , J 7 h-> M 3 , Js ^ Mi] 



5 



x> 



-• • 



IJi> 



|0> Pi1 Pi2 Pi3 Pi4 



Figure 2: A quantum circuit to prepare the job state vectors with m = 4. In this case | x) is a set of 
two control qubits. Each control qubit is set to l/\/2(| 0)+ | 1)). The target qubits will be prepared 
in a superposition of all possible job-machine states. 
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Figure 3: A quantum circuit to prepare | Ji). The execution times of J\ are 1, 3, 7, and 15 units of 
time respectively. | x) are control qubits which control the selection of the four job-machine states, 
e ii 1 < j < 4. The first two qubits of | Ji) record the machine index, and the remaining four qubits 
record the time if Ji is assigned to the corresponding machine. 



Schedule Si corresponds to the following job state vectors: 

| Ji) =| 000001) | J 2 ) =| 010001) | J 3 ) =| 000H0) | Ji) =| 110100) 

I Jh) =| 100011) | J 6 ) =| 010111) | J 7 ) =| 100011) | J 8 ) =| 000001) 

The schedule vector for schedule Si is: 

| Si) =| 000001)® | 010001)® | ooono)® | 110100) 

<g> | 100011)® | 010111)® | 100011)® | oooooi) 
The completion times on all machines are: 

C Al (Mi) = l + 6 + l = 8, C Al (M 2 ) = l + 7 = 8, C Al (M 3 ) = 3 + 3 = 6, C Al (M 4 ) = 4. 

The makespan of schedule Si is equal to the largest completion time over all machines: 



n 

(ii) Schedule S 2 : 



At 

max 
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[Jx^M.2, Jz^Mx, J 3 ^M 3 , Ja^M 4 , J 5 ^M 3 , J 6 ^M 2 , Ji ^ Mi, J% h-> M 2 ] 



The schedule vector is: 

I S 2 ) =| 010011)® | 000010)® | 100101)® | 110100) 

<8> | 100011)® | 010111)® | 000101)® | 011010)) 

The schedule vector can also be in a superposition of some basic states, for example, the schedule 
vector could be in an equal superposition of the two schedules, Si and S 2 : 

| S) = 1/V2(| Si) + | S 2 )) 

= l/\/2(| 000001)® | 010001)® | 000110)® I 110100) 

® I 100011)® I 010111)® I 100011)® I 010001) 
+ I 010011)® I 000010)® I 100101)® I 110100) 

® I 100011)® I 010111)® I oooioi)® I 011010)). 

4 The Running Time of Machine A4j under Schedule Sk 

Computing the total running time of machine Mj under schedule Sk requires summation of the 
running times of all jobs assigned to Mj. A quantum adder similar to a classic adder, e.g. | 5)+ | 
6) =| 11), has: 

• n inputs oi, a 2 , • • ■ , a n , each one is a register of q qubits, 

• one carry in c, and 

• one output S = Y17=l Q-i + c, a register of q + n qubits. 

To obtain the total running time of machine Mj under schedule Sk we add the execution times, 
(in our examples those qubits of a job-machine vector which are not underlined) for all jobs assigned 
to Mj and create a running time vector for machine Mj, \ Tj). The running time vector for schedule 
Sk is the tensor product of the execution time on individual machines under schedule Sk- 

\t s *)=\t*><)®\t£*)---®\t^). 

For example, for machine Mi under schedule Si (see Section 0): 

| Tf 1 ) : | 0001)+ | 0000)+ | 0110)+ | 0000)+ | 0000)+ | 0000)+ | 0000)+ | 0001), 

or 

| Tf 1 ) : | 0001000). 

We want to construct a quantum circuit to sum the execution time of jobs assigned to each 
machine Mj . We use the index qubits as control qubits to control the summation of each index, 
or each machine; the index qubits are entangled with the target qubits which give the running time 
on that machine. As the input qubits | S) are prepared in the equal superposition of all possible 
schedules, this operation will prepare the running time of machines under all possible schedules. 

In Figure^ T\, T 2 , . . . Tm represent the execution time on machines Mi,M 2 , ...Mm respectively. 
When the index qubits are not "active" , the respective input for the summation circuit will be zero. 
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Figure 4: A circuit to compute the sum of the execution time of jobs assigned to each machine. 
Jl, Ji . . . Jjv are the job vectors prepared according to the scheme discussed in Section |3 and 
Ti,T2, . . .Tm represent the execution time on machines A4i,M 2 , ...M.m respectively. 



For example, if the index of e| =| 010001) is 01, rather than 00, then the index qubits of e\ for | Ti) 
are not active. 

How to implement arithmetic on a quantum computer has been address by many authors. Many 
detailed quantum circuits and discussions can be found in |71 151 IT3]. For the example discussed in 
Section |21 after the summation operation, the system will be in the state: 



| ST) = l/V2{\S l T 1 )+\S 2 T 2 )) 

| ST) = l/\/2(| 000001)® | 010001)® | 000110)® | 110100)® | 100011)® | 010111) 

® | 100011)® | 010001)® I 0001000)® I 0001000)® I 0000110)® I 0000100) 

+ I 010011)® I 000010)® I 100101)® I 110100)® I 100011)® I 010111) 

® 1 000101)® I 011010)® I 0000111)® I 0010100)® I 0001000)® I 0000100)) 



5 Determination of the Makespan 



Now we have the system in state | ST), of which the last M(n + q) qubits provide the running 
time of all M machines under all schedules. The makespan of a schedule is equal to the maximum 
running time among the M machines under that schedule. We want to construct a Max circuit, as 
shown in Figure UJ to compute the makespan of each schedule. The quantum circuit computes the 
maximum over an input set, e.g., Max(| 5), | 7), | 4)) =| 7). The input to this circuit is a set of n + q 
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qubits. The output | Max), has also n + q qubits. Implementing such arithmetic on a quantum 
computer has been addressed by many researchers [71 151 114) and we omit the detail circuit here. 



Ti 
T 2 



Tm 
|0> 



✓max 



Figure 5: A quantum circuit to compute the makespan, the maximum running time among all 
machines. The output of this circuit is an entangled state. 

The output of a quantum circuit as in Figure |S] is an entangled state, rather than the tensor 
product | Ti)(g) | T2) ® • • • ® | Tm)® \ C ma x}- Recall that the | T) was prepared in an equal 
superposition of running times of the machines under all possible schedules. Thus, such a Max 
operation prepares the C max register in an equal superposition of the makespans of all possible 
schedules. 

Up to now, we discussed the implementation of the quantum circuits presented in Figure ^ 
During these operations, we successfully entangle the job vectors with the Sum and the Makespan 
vectors. These vectors are prepared in the equal superposition of all possible schedules, which can 
be written as: 



' " x 

STC max ) — y ^ SiTiCn 



1 



where, 1 < i < 2 m2 indexes different schedules, and we ignore the coefficients. 
For the simple example in Section the system will be in the state : 



STC max ) = 1/V2(| 50 I T x ) I C maxl )+ I S 2 ) I r 2 ) 1 c max2 )) 

= V\/2(| 0Q0001)® 1 010001)® 1 000110)® 1 110100)® 1 100011)® 1 010111)® 1 100011) 
® 1 010001)® 1 0001000)® 1 0001000)® 1 0000110)® 1 0000100)® 1 0001000) 

+ 1 010011)® 1 000010)® 1 100101)® 1 110100)® 1 100011)® 1 010111)® 1 000101) 

® 1 011010)® 1 0000111)® 1 0010100)® 1 0001000)® 1 0000100)® 1 0010100)) 



6 Searching for a Schedule with a Given Makespan 

In our approach, a generalized version of Grover's search algorithm allows us to find a schedule with 
a given makespan. The basic ideas of Grover's quantum search algorithm are discussed next. 
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Consider a search space T searc h = {E x } consisting of N — 2™ elements. Each element E x , 1 < 
x < 2™, is uniquely identified by a binary n-tuple x, called the index of the element. We assume 
that M < N elements satisfy the requirements of a query and we wish to identify one of them. 



n 

qubits 

oracle 
workspace 



oracle 



|x> ^(-1) f(x) |x> 




V 

U = 2 |0><0| - / 



Figure 6: A quantum circuit for the Grover's iteration 

The classic approach is to repeatedly select an element Ej, decide if the element is a solution 
to the query, and if so, terminate the search. If there is a single solution (M = 1) then a classical 
exhaustive search algorithm requires 0(2") iterations. 

For the Grover quantum searching algorithm, we apply a Walsh-Hadamard transform to create 
an equal superposition state which includes all elements of the search space: 

N-l 

i </>> = -7= y i x). 

v x— 

Then we perform Grover's iterations. The circuit for this algorithm is presented in Figure 

An oracle examines an index/label, x, and decides if it matches the search argument or not. To 
abstract this process we consider a function f(x) with < x < 2™ — 1 such that 

, . JO if a; is not a solution 
^ ' 1 1 if a; is a solution. 

An oracle qubit, \ q), initially set to | 0) is reset to | 1) when the oracle recognizes a solution to the 
search problem we pose. The black box oracle O performs the following transformation 

0\x)\q) =\x)\q®f(x)). 

The oracle qubit can be initially in the state | q) = {l/^/2){\ 0)— | 1)). 
Thus this transformation can be rewritten as 

0\x){\ 0)- | 1»/V2= | x) (| 0)- | l))A/2. 

The state of the oracle qubit does not change and can be omitted from the description of the quantum 
search algorithm 

| x) i — ^ (-l/( x > | x). 
Let U be the following transformation: 

f/ = 2|0)(0|-J. 
Then a conditional phase shift in Figure applied to the system is: 

S p = H® n UH® n = H® n {2 | 0)(0 | -I)H® n =2\^){ij}\ -I. 
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A Grover's iteration consists of O, the transformation performed by the oracle followed by a condi- 
tional phase shift: 



G = S p O = (2 | I -1)0 
Thus, the quantum search algorithm could be written as: 

N-l -, N-l 



GR ^f E I x ) I ?> = K2 I 1>)(1> I ~W4= E I x K\ °>" I x )(| 0)- | 1))/V2 



X — a;— u 

When after i? = 0{\J~^) iterations, we measure the first n qubits and obtain xq, the solution to 
the search problem. 

Amplitude Amplification represents a generalization of the Grover's quantum search idea |2*llll|. 
Let A be a unitary operator in a Hilbert space, Hn with an orthonormal basis | 0), | 1), . . . | N — 1); 
the only condition imposed on A is to be invertible, thus A must not involve any measurements. 

If X : {0, 1, ... AT — 1} i—> {0, 1} is a Boolean function we say that the basis state | x) is a "Good" 
state if x{ x ) — 1 an d | x) is a "Bad" state" if x{ x ) = 0. The central piece of the amplitude 
amplification is an operator Q defined as: 

Q=Q(A,x,4>,ip) = -AS (<t>)A- 1 S x (<p). 

with cf> and tp two angles such that < 4>, ip < 7r and S x an operator which conditionally changes 
the amplitudes of "Good" states: 

e iv | x) if x(x) = 1 
\x) if x(x)=0. 

Similarly, So amplifies the amplitude by a factor e 4 * if the state is not | 0) . 

Let a denote the probability of finding a "Good" element x; amplitude amplification allows to 
find a "Good" x after an expected number of applications of A and of the inverse of A; the number 
of iterations is proportional to l/y/a. We also define the angle 9 such that: 

sin 2 (<9) = a. 

Grover's algorithm is a particular instance of amplitude amplification when the oracle implements 
the Boolean function / = x, and the transformation A is the Walsh-Hadamard transform W = H® n 
on n qubits. 

This iteration carried out by transformation Q can be regarded as a rotation in the two- 
dimensional space spanned by the the state of a uniform superposition of non-solutions and the 
state consisting of a uniform superposition of solutions to the search problem. The initial state may 
be expressed as: 



| ^o) = -A I Good) + VI - a | Bad) 
Figure [7| presents the effect of the transformation Q = —ASqA~ 1 S x as: 

• the oracle operation S x performs a reflection about the vector | Good) . 

S x | (x(x) = l) S x \x)=-\x) (x(x)=0) 

• ASqA^ 1 performs a reflection about the initial state | ipo) 

So | 0) HO} S \x) = -\x) (x?0) 
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|Good> 



|Good> 



|Good> 




Figure 7: The search operator Q performs a rotation toward | Good) states by 29 radians. 8 is the 
angle between the initial state | ipo) and the | Bad) in the two-dimensional space spanned by the 
| Good) and | Bad) states and sin 2 (6*) = a with a the probability of finding a "Good" element x. (a) 
The current state | ip) and the initial state \ ipo)- (b) The oracle operation S x performs a reflection 
of the current state | ip) about the vector | Good), (c) ASqA^ 1 performs a reflection about the 
initial state | ipo) 



• Q = — ASqA 1 S x performs a rotation toward | Good) vector by 29 radians, where sin 2 9 = a 

Each iteration Q will rotate the system state by 28 radians toward the solutions of the searching 
problem. Thus after m iterations, the measurement on the final state Q m \ ipo) will produce a 
"Good" state with probability equal to sin 2 ((2m + 1)9). The amplitude amplification algorithm 
could find a "Good" solution in 0(-^) 0^5- Each iteration involves the application of A and 

A- 1 . 

Now let us return to our scheduling problem and recall that: 

• We prepare each job vector | Ji) in a superposition state which includes the running times on 
all machines. The first m qubits of a job vector are used for the index of the machine and 
remaining q qubits are used for the running time of that job on the machine. 

• We summarize the execution time of all jobs according to their machine indexes, which produces 
all schedules and gives the running time of each machine Tj under these schedules. The system 
is prepared in an entangled state: 

N M 

|5T) = -= Y, (®l-**>®l r i*>) <t = 2™ 2 ", 

each schedules k i— 1 j — 1 

a superposition of job vectors and running time vectors of all possible schedules. 

• We obtain the maximum running time among all machines using the Max quantum circuit and 
prepare the system in state: 



stcZ2) = 4= Yl (® I J ^ ® I T ^ I C ™°* fe » 



N M 

m2" 



each schedules k i— 1 j — 1 



As we can see, our algorithm successfully prepares the system in an equal superposition of all 
2 m2 possible schedules. We define this whole preparation process as Z. This Z transformation 
does not carry out a measurement of the system at any time. Therefore, there exists an inverse 
transformation operation Z~ x . We can use the amplitude amplification algorithm to search the 
schedule with the makespan Dk = /i. If we find such a makespan, the job vectors will be projected 
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as well. These projections will give us the actual mapping of jobs to machines corresponding to the 
schedule with the given makespan. 

The searching process consists of the following steps: 

• Apply the Z transformation on | 0) to prepare the system in an equal superposition of all 2 m2 
possible schedules, | tp). 

• Repeat the following steps 0(y/a) times: 

Apply Q on | tp), in which Q = —ZSqZ~ 1 S x 

• Measure the resulting state. 

• Return the result; the job vectors give the detailed schedule. 

The oracle in our searching algorithm exhibits some difference with the oracle in Grover's al- 
gorithm which checks all qubits to reverse the solution(s) of the searching problem. In our case, 
the oracle only checks a subset, C max qubits. It is easy to implement such an oracle using an idea 
similar to the one in |17| . 

The algorithm presented in this paper can be optimized; to avoid over rotation we could use the 
fix-point quantum search |1U| . 



7 Scheduling Problems with a Quantum Counterpart 

Many other scheduling problems can be reformulated to take advantage of quantum search and 
exhibit a square root speedup versus their classical counterparts. Such problems require that a 
synthetic measure of performance, /i be within a given range, fi m in < A* < Umax- 
The generic quantum algorithm proceeds as follows: 

1. Devise an encoding scheme for the information required to compute for a given schedule 
Si. 

2. Design an algorithm to compute the synthetic measure of performance, p,s i - 

3. Construct | Q) with Q = (/^s i ,S'i) in a superposition state for all possible schedules. 

4. Design the quantum circuits for the specific encoding scheme and for the algorithms. Given a 
set of values {q%, 92, ■ ■ ■ In} the circuit should be able to compute a function fig t = f(q±, (72, • ■ ■ ?n)- 
For example, we may wish to compute the maximum, the minimum, or an average value. 

5. Design an oracle to identify a specific value of the function fig i . 

6. Use the quantum search to find if there is a value ns i = fi m in- If so, determine the correspond- 
ing schedule Si. Continue this process until fj,s i — fAmax. 

7. If no schedule can be found then report failure, otherwise provide the list of all schedules Si 
and the corresponding measures of performance, fj,s i ■ 

Consider for example the -R||5Zi= 1 scheduling problem when the goal is to optimize the 
average completion time in unrelated parallel machine environment. It has the similar encoding 
process as the R\\C m ax problem. During the process we construct all schedules, we also summarize 
the completion time of different jobs using some simple arithmetic circuits, followed by Grover-type 
search. Other scheduling problems such as minimizing the average waiting time, could also take 
advantage of quantum search. 

Oftentimes, we have to search for schedules that optimize the largest subset of a set of synthetic 
measures of performance, fi, v, tt, p, For example, we could have multiple synthetic performance 
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indicators related to: timing, resource utilization, cost, and quality of service. In this case we would 
run repeatedly the scheduling algorithm for each performance measure and search for a solution in 
the given range for each measure. Once we have conducted all individual searches we determine 
the intersection of all schedules that satisfy all conditions; if the set is empty we remove individual 
conditions one by one until we find the a non-empty set. 

Scheduling is also intimately related to planning when we have a complex goal and the sequence 
of actions to reach each goal have to be determined. The scenario described above is repeated for 
each plan thus the square root speedup of algorithms based upon quantum search becomes even 
more important. 



8 Summary 

When a deadline is imposed, or when we wish to find a schedule with a given range of possible average 
completion time we discover that a full range of scheduling problems have a quantum counterpart 
which can take advantage of Grover's search. 

Many scheduling problems, resource allocations, and path-finding problems, share some common 
properties with the R\\C max problem discussed in this paper: a well-defined initial state, a well- 
defined desired state or a range of desired states, many paths to reach the desired state, and well- 
defined attributes of an optimal path. 

The quantum algorithmic solution to such problems requires the following steps: 

• Prepare the initial state in an equal superposition of all possible choices. 

• Use some reversible quantum arithmetic to compute the specialized property (makespan in 
our case) needed. 

• Construct the necessary oracle circuit. 

• Use Grover-type algorithms to search for the desired solution. 

The solution we propose based upon Grover's algorithm is not universally applicable. Problem 
requiring specific optimization criteria may require quantum circuits that cannot be easily imple- 
mented with reversible quantum arithmetic. Moreover, Grover-type algorithms lead to a square-root 
improvement over the exhausting search, while many good classical algorithms may have better per- 
formance for some special problems. 
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